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97.16  Probability  analysis  for  rolls  of  a  square  cuboidal  die 

Consider  a  six-sided  die  in  which  two  opposite  faces  are  squares  of  sides  y 
and  the  other  four  faces  are  rectangles  of  sides  x  and  y,  as  sketched  in  Figure  1 . 
One  can  think  of  it  as  an  ordinary  cubical  die  that  has  been  stretched  or 
compressed  along  one  axis.  There  is  some  probability  p  that  one  of  the  two  square 
faces  will  show  uppermost  after  such  a  die  is  randomly  tossed,  and  a  probability 
I  -p  that  one  of  the  four  rectangular  faces  will  show.  For  example,  in  the  limit 
that  y  becomes  equal  to  x  so  that  the  die  is  cubical,  then  the  probability  that  the 
original  square  faces  come  up  is  p  =  5.  Over  30  years  ago,  the  question  was 
asked  f  1]  what  is  p  if  y  is  not  equal  to  x?  The  author  of  that  article  cut  15  such  dice 
out  of  a  steel  bar  and  had  his  students  roll  each  of  them  N  times  to  experimentally 
determine  p.  His  results  are  listed  in  Table  1.  A  simple  explanation  was 
subsequently  hypothesized  in  which  p  is  proportional  to  the  solid  angle  subtended 
by  a  square  face  [2],  but  it  does  not  match  the  data  well.  In  the  present  paper,  an 
alternative  simple  model  is  proposed  that  better  fits  the  measurements. 
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FIGURE  1  FIGURE  2 


The  model  is  developed  in  three  steps.  First  suppose  the  die  lands 
vertically  without  initial  rotation  on  any  one  of  its  eight  edges  of  length  y  in 
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Figure  1 .  Let  the  bottom  square  face  make  angle  0  with  respect  to  the  floor 
at  the  instant  of  landing,  as  sketched  in  Figure  2.  The  die  will  subsequently 
fall  over  either  onto  that  bottom  face  or  onto  the  right-hand  rectangular  face, 
depending  on  whether  the  centre  of  mass  of  the  cuboid  is  to  the  left  or  the 
right  of  the  edge  that  initially  hits  the  floor  [3],  In  other  words,  it  depends  on 
whether  a  diagonal  drawn  across  the  front  face  of  the  die,  as  shown  in 
Figure  2,  is  inclined  leftward  or  rightward  of  vertical.  That  diagonal  makes 
angle  xp  =  tarT1  (x/y)  relative  to  the  bottom  face  and  it  will  be  exactly 
vertical  when  the  die  makes  angle  60  =  tt/2  -  (p  relative  to  the  floor.  The 
die  will  fall  onto  the  square  face  if  0  <  6  <  60.  Assuming  that  all  landing 
angles  0  between  0  and  jr/2  are  equally  likely,  the  probability  pm  that  the 
die  ends  with  the  square  face  up  in  this  two-dimensional  case  without  initial 
rotation  is  therefore 
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Intuitively  it  makes  sense  that  the  probability  only  depends  on  the  ratio  of  x 
and  y,  and  not  on  those  two  lengths  independently. 


The  second  step  is  to  extend  this  result  to  find  the  probability  pjD  of 
ending  with  a  square  face  up  if  the  die  lands  vertically  without  initial 
rotation  on  any  of  its  twelve  edges  as  the  cuboid  is  reoriented  in  three 
dimensions.  The  much  more  likely  case  that  it  lands  on  a  corner  can  be 
thought  of  as  a  superposition  of  landing  on  the  two  nearest  lower  edges  of 
the  die  and  is  therefore  included  in  py>.  (If  the  corner  of  a  large  box  is  held 
against  the  floor  and  the  box  is  then  released,  separate  rotations  about  each 
of  the  lower  edges  are  observed  as  the  box  settles.)  If  the  die  in  Figure  1 
lands  on  one  of  the  four  edges  of  length  x,  it  will  end  with  a  rectangular  face 
up.  So  pio  can  be  calculated  by  weighting  pin  by  the  fractional  probability 
8y/  (8y  +  4x)  that  the  die  (which  spins  about  a  random  axis  while  in  the  air) 
will  land  on  an  edge  of  length  y  instead  of  one  of  length  x, 

Pv>  =  o-  8>  *P2D  =  „  T1  -  “  tan-1  (-)l .  (1) 


8y  +  4x  ~  2  +  x/y  [_  tt  \y) J 

This  equation  correctly  predicts  py,  =  1  for  a  cubical  die  with  x  =  y. 


The  final  step  is  to  incorporate  the  fact  that  when  a  real  die  is  thrown,  it 
lands  on  the  table  with  substantial  initial  horizontal  translational  and  angular 
velocities.  Consequently  it  tumbles  along  the  table  for  a  while.  The  die  slows 
down  each  time  a  face  inelastically  hits  the  table  [4],  It  can  be  supposed  that 
the  loss  is  proportional  to  the  area  A  of  the  colliding  face  (and  hence  the  die  is 
even  more  likely  to  end  its  motion  on  a  larger  area  face)  for  two  reasons.  First, 
a  greater  number  of  elastic  modes  of  vibration  can  be  excited  over  a  bigger 
surface.  (As  proof,  a  block  of  wood  with  an  aspect  ratio  x/y  of  about  10  was 
dropped  onto  a  concrete  floor  from  a  height  of  about  1  m.  It  bounced 
noticeably  higher  when  dropped  on  its  end  compared  to  on  its  side.)  Secondly, 
the  larger  faces  of  the  cuboid  will  hit  the  table  with  higher  impact  speed  (and 
hence  loss)  than  the  smaller  faces  as  a  die  rolls  around  any  edge  y  because  the 
die  has  lower  centre  of  mass  in  the  former  orientations.  To  model  this  rolling 
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loss,  the  factors  of  x/y  in  (1)  can  thus  be  weighted  by  some  positive  power  of 
the  ratio  of  the  rectangular  and  square  areas,  xyly 2,  to  give  a  final  prediction  of 
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where  n  >  1.  The  exponent  n  effectively  makes  a  thin  die  (for  which  x  <  y) 
thinner  and  a  thick  die  (x  >  y)  thicker.  For  example,  if  n  =  3  then  a  rolled  die 
whose  rectangular  edge  x  is  half  the  length  of  its  square  edge  y  will  behave  like 
a  die  with  x  =  y  1 8  that  lands  vertically  without  initial  rotation.  This  example 
suggests  that  the  value  of  n  may  depend  on  the  height  and  flick  of  the  wrist 
used  to  throw  the  die.  One  might  also  expect  n  to  depend  on  the  coefficient  of 
restitution  (that  is,  on  the  hardness  of  the  die  and  table),  on  any  slight  beveling 
of  the  edges  of  the  die,  and  on  similar  factors.  In  support  of  this  expectation,  a 
second  set  of  data  for  square  cuboidal  dice  made  of  polyvinylchloride,  rather 
than  steel,  has  been  collected  by  Heilbronner  [5]  as  listed  in  Table  2;  they  have 
somewhat  different  values  of  p,  presumably  because  of  their  different 
compositions.  It  is  evidently  necessary  that  any  general  model  such  as  (2)  have 
some  empirical  parameter  in  it.  For  example,  another  way  to  account  for  the 
reduced  likelihood  of  the  die  ending  up  in  an  orientation  with  high  centre  of 
mass  is  to  weight  the  probability  by  a  Boltzmann  exponential  exp  (-/?£),  where 
E  is  the  gravitational  potential  energy  of  the  centre  of  mass  of  the  die  and  /?  is 
taken  to  be  a  fit  parameter  [5,  6,  7]. 
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The  value  of  the  parameter  n  is  found  by  fitting  (2)  to  the  data.  The  results 
are  graphed  in  Figure  3  for  two  values  of  n  that  span  the  range  of  variation  of 
both  Budden's  and  Heilbronner's  experimental  points.  The  conclusion  is  that 
n  =  3.0  ±  0.5  for  cuboidal  dice  rolled  in  an  ordinary  manner.  (Consequently 
even  a  slight  shaving  of  one  face  of  a  die  would  significantly  alter  its  fairness.) 
This  same  model  could  be  applied  to  the  related  problem  of  the  probability 
1  -  p  for  a  tossed  coin  of  thickness  x  and  diameter  y  to  land  on  edge  [2,  4,  8,  9] 
if  one  approximates  a  coin  as  behaving  like  the  smallest  cuboid  that  would 
enclose  it.  In  particular,  (2)  predicts  the  probability  is  between  0.12%  and 
0.83%  for  a  British  one-pound  coin  (of  3.15  mm  thickness  and  22.5  mm 
diameter)  if  n  is  between  2.5  and  3.5;  this  range  brackets  the  experimental 
observation  of  6  edge  landings  out  of  1000  tosses  [8]. 
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- n  =  2.5  - n  =  3.5  •  Budden  experiment  ■  Heilbronner  experiment 
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FIGURE  3 
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